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Abstract. 

A detailed analysis of the radiative and Auger de-excitation channels of K-shell vacancy states in Fe n-Fe ix has been carried out. 
Level energies, wavelengths, A-values, Auger rates and fluorescence yields have been calculated for the lowest fine-structure 
levels populated by photoionization of the ground state of the parent ion. Different branching ratios, namely Ka 2 /Ka u K0/K cr, 
KLM/KLL, KMM/KLL, and the total K-shell fluorescence yields, co K , obtained in the present work have been compared with 
other theoretical data and solid-state measurements, finding good general agreement with the latter. The Kcr^/Kc*! ratio is found 
to be sensitive to the excitation mechanism. From these comparisons it has been possible to estimate an accuracy of -10% for 
the present transition probabilities. 
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1. Introduction 

The iron K lines appear in a relatively unconfused spectral re- 
gion and have a high diagnostic potential. The study of these 
lines has been encouraged by the quality spectra emerging from 
Chandra and by the higher resolution expected from Astro-E 
and Constellation^ . In addition there is a shortage of accurate 
and complete level-to-level atomic data sets for the K- vacancy 
states of the Fe isonuclear sequence, in particular for the lowly 
ionized species. This undermines line identification and real- 
istic spectral modeling. We are currently remedying this situa- 
tion by systematic calculations using suites of codes developed 
in the field of computational atomic physics. Publicly available 
packages have been chosen rather than in-house developments. 
In this context, complete data sets for the n = 2 K- vacancy 
states of the first row, namely Fe xvm-Fe xxv, have been re- 
ported earlier by Bautista et al. (2003) and Palmeri et al. (2003), 
to be referred to hereafter as Paper I and Paper II, and for the 
second row (Fe x-Fe xvii) by Mendoza et al. (2003), to be re- 
ferred to as Paper III. 

The K lines from Fe species with electron occupancies 
N > \1 have been studied very little. Jacobs & Rozsnyai 
(1986) have computed fluorescence probabilities in a frozen- 
core approximation for vacancies among the nt subshells of 
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the Fe isonuclear sequence. Kaastra & Mewe (1993) have cal- 
culated the inner-shell decay processes for all the ions from Be 
to Zn by scaling published Auger and radiative rates in neu- 
trals. Both of these studies ignore multiplets and fine-structure. 
Otherwise, the bulk of the data in the literature is devoted to 
solid-state iron. Regarding wavelengths and line intensity ra- 
tios, numerous references are listed in Holzer et al. (1997) who 
measure the Ka\ t2 and K£i ,3 emission lines of the 3d transi- 
tion metals using a high-precision single-crystal spectrometer. 
Fewer publications are available on the experimental K Auger 
spectra: Kovalfk et al. (1987) have derived the KLM/KLL and 
KMM/KLL ratios from the Auger electron spectrum of iron 
produced by 57 Co decay, and Ndmethy et al. (1996) have mea- 
sured the KLL and KLM spectra of the 3d transition metals but 
have not determined the KLM/KLL ratio. Concerning K-shell 
fluorescence yields, measurements covering the period 1978- 
93 have been reviewed by Hubbell et al. ( 1 994) following major 
compilations by Bambynek et al. (1972) and Krause (1979). 

The present work is a detailed analysis of the radiative and 
Auger de-excitation channels of the K-shell vacancy states in 
the third-row species Fe ii-Fe ix. Energy levels, wavelengths, 
A-values, Auger rates and fluorescence yields have been com- 
puted for the lowest fine-structure levels in configurations ob- 
tained by removing a Is electron from the ground configura- 
tion of the parent ion. In Section 2 the numerical method is 
briefly described. Section 3 outlines the decay trees and se- 
lection rules. Results and discussions are given in Section 4, 
followed by a summary and conclusions (Section 5). All the 
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atomic data calculated in this work are available in the elec- 
tronic Tables 3^. 


2. Numerical method 


In Paper III we conclude that the importance of core-relaxation 
effects increases with electron occupancy, so calculations of 
third-row iron ions require a computational platform which is 
well suited to treating these effects. For this reason the cal- 
culations reported here are carried out using the hfr pack- 
age by Cowan (1981) although autostructure (Badnell 1986, 
1997) is also heavily used for comparison purposes. In hfr 
an orbital basis is obtained for each electronic configuration 
by solving the Hartree-Fock equations for the spherically av- 
eraged atom. The equations are derived from the application 
of the variational principle to the configuration average en- 
ergy and include relativistic corrections, namely the Blume- 
Watson spin-orbit, mass-velocity and the one-body Darwin 
terms. The eigenvalues and eigenstates thus obtained are used 
to compute the wavelength and A-value for each possible tran- 
sition. Autoionization rates are calculated in a perturbation the- 
ory scheme where the radial functions of the initial and final 
states are optimized separately, and configuration interaction 
is accounted for only in the autoionizing state. Configuration 
interaction is taken into account among the following config- 
urations: (3d + 4s) w , [3p](3d + 4s) Af+1 , [2p](3d + 4s) M+1 and 
[ls](3d + 4s) M+1 where [i nt ] stands for a hole in the nt sub- 
shell and {nt + n't') M represents all possible distributions of M 
electrons among the nt and n't' shells, M ranging from M - 0 
in Fe ix to Af = 7 in Fe n. 

Given the complexity of [2p][3p] double- vacancy channels 
in ions with an open 3d shell, the level-to-level computation of 
the Auger rates with hfr and autostructure proved to be in- 
tractable. However, average values can be used for all the lev- 
els taking advantage of the near constancy of the total Auger 
widths in ions for which the KLL channels reduce the outer 
configuration to a spectator (see Paper III). Therefore, we em- 
ploy the formula given in Palmeri et al. (2001) for the single- 
configuration average (SCA) Auger decay rate 


sca = L(27 t + lMa(0 

‘ L(2 Ji + 1) 


—A a {ntriC -4 n"i"et"‘) 
g 


( 1 ) 


where the sum runs over all the levels of the autoionizing con- 
figuration, s/g is a statistical factor given in Eqs. (15-16) of 
Palmeri et al. (2001) that contains the dependence on the active 
shell ( nt , rit\ n"t") occupancy, and A a (ntn't' -> ri’t”et f ") is 
the two-electron autoionization rate which is a function of the 
radial integrals and for which the complete expression is given 
in Eq. (1 1) of Palmeri et al. (2001). The calculated Auger rates 
are expected be as accurate as those obtained in a level-to-level 
single-configuration approach. 


3. Photoionization selection rules and decay trees 


We have focused our calculations on the Fe K- vacancy states 
populated by photoionization of the ground state of the parent 
ion 


// ^Lj+y 


( 2 ) 


where // is the outer configuration of the ground state, i.e. 
3d 6 4s 2 in Fe i, 3d 6 4s in Fe n and 3d jV ( N = 1 - 6) in Fe in— 
Fe vm, [Is] denoting a K hole. In Eq. (2), K-shell photoion- 
ization leads to a p wave with selection rules specified by Rau 
(1976): L' = L, S' - S = ±1/2 and J' - J = ±1/2, ±3/2. 
Therefore, only few states are expected to be populated. 

The radiative and Auger decay manifolds of a K- vacancy 
configuration [Is ]/i can be outlined as follows: 


Radiative channels 

Kfi 

[ls]/i — > [3p \ji + yp 
^ [2p]/i + y a 
Auger channels 


[Is 


KMM(p) 


KMM(s) 


KLL 


(3) 
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[3s] 2 // + e~ J 

(6) 

[2s]// -1 + e~ \ 
[2p]// -1 +e-j 

(7) 

[2s] [3p] n + e~ 1 
[2p] [3p]// + e~ 
[2s] [3s] /i + e~ 
[2p] [3s] fx + e~ 
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where the negative exponent of /i stands for the number of elec- 
trons that have been extracted from the outer subshells. In the 
radiative channels, forbidden and two-electron transitions have 
been excluded as it has been confirmed by calculation that 
they display very small transition probabilities (logA r < 11). 
Therefore, the two main photo-decay pathways are the result 
of the 2p —> Is and 3p -» Is single electron jumps that give 
rise respectively to the Kor (~ A1.94) and Kp (~ /11.75) arrays. 

It has also been numerically verified that the participator 
Auger channels (KMM(p) and KLM(p)) contribute less than 
one percent to the total Auger widths, and hence, they have not 
been taken into account. Of primary interest in the present work 
is the branching ratios of the Kor, K/3, KMM, KLM and KLL 
channels and their variations with electrons occupancy N. 


4. Results and discussion 

Centroid wavelengths for the K a and Kfi unresolved transition 
arrays (UTAs) computed with hfr are presented in Table 1 , in- 
cluding also a comparison with experiment. For the Kct\ and 
K a 2 lines, it can be seen that the agreement with the solid-state 
measurements by Holzer et al. (1997) (~ 0.5 mA) is somewhat 
better than that with the EBIT results for Fe x (Decaux et al. 
1995) of within 2 mA, and the slight blueshift with increas- 
ing N predicted by hfr is consistent with experiment. A small 
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Fig. 1. Comparison of centroid wavelengths for (a) K a and (b) Kfi 
UTAs in Fe ions with 17 < N < 25. Filled circles: this work. Filled 
triangles: Kaastra & Mewe (1993). 


Table 1. Comparison of centroid wavelengths (A) for the K a and 
Kfi unresolved transition arrays in Fe ions with 17 < N < 25. The 
hfr wavelengths have been weighted with the fluorescence yields. 
Experimental wavelengths for Fe x (N = 17) are from Decaux et al. 
(1995) and for solid-state iron from HSlzer et al. (1997). 



HFR 

Expt 

N 

Ka 

W 

Ka K/? 

17 

1.9369 

1.7488 

1.9388(5) 


1.9407 

1.7493 

1.9413(5) 

18 

1.9366 

1.9402 

1.7511 


19 

1.9364 

1.9403 

1.7527 


20 

1.9362 

1.9405 

1.7540 


21 

1.9360 

1.9404 

1.7557 


22 

1.9357 

1.9400 

1.7566 


23 

1.9355 

1.9397 

1.7565 


24 

1.9356 

1.9398 

1.7565 


25 

1.9356 

1.9398 

1.7566 

1.936041(3) 1.756604(4) 

1.939973(3) 



Fig. 2. Comparison of K/?/K a intensity ratios for Fe ions with 18 < 
N < 25. Filled circles: hfr, this work. Dotted upright triangles: au- 
tostructure, this work. Dotted diamonds: mcdf-sal, this work. Filled 
upright triangles: Kaastra & Mewe (1993). Filled squares: Jacobs & 
Rozsnyai (1986). Filled diamonds: Scofield (1974). Filled inverted 
triangles: Jankowski & Polasik (1989). Circles: Perujo et al. (1987). 
Upright triangles: Holzer et al. (1997). Squares: Rao et al. (1986). 
Diamonds: Berdnyi et al. (1978). Inverted triangles: Salem et al. 
(1972). Dotted circles: Sli vinsky & Ebert (1972). Dotted squares: 
Hansen et al. (1970). 


Table 2. Comparison of hfr Ka 2 /Ka\ ratios for Fe ions ( 1 8 < N < 25) 
with experiment and previous theoretical estimates. 


N 

Ion 

HFRl a 

HFR2 b 

Expt. 

Theory 

18 

Feix 

0.7000 

0.5912 



19 

Fe vm 

0.5780 

0.6015 



20 

Fe vn 

0.4541 

0.5821 



21 

Fe vi 

0.5545 

0.5538 



22 

Fe v 

0.5852 

0.5346 



23 

Feiv 

0.4426 

0.5143 



24 

Fe m 

0.4438 

0.5032 



25 

Fen 

0.4476 

0.4969 

0.5 1 (2) c 
0.507(10) d 
0.506 e 
0.4998 f 

0.5107* 


a Includes only the decay lines from the K-vacancy levels listed in 
Table 3 

b Includes all the decay lines from the [ls](3d +4s) w+1 complex 
c Holzer etal. (1997) 
d McGrary etal. (1971) 
e Salem &Wimmer (1970) 
f Williams (1933) 

1 Scofield (1974) 


redshift with N is also found for the K/? array. On the other 
hand, as shown in Figure 1 , the present findings contrast with 
the steeper redshift for both UTAs obtained from the data by 
Kaastra & Mewe (1993). 

In Table 2, Ka 2 /Kct\ intensity ratios are tabulated for differ- 
ent ionization stages. Two cases have been considered: HFR1, 
only the decay lines from the K-vacancy levels populated by 
the photoionization of the respective ground state are included; 
and HFR2, all the transitions from the levels belonging to the 
[ls](3d + 4s) w+1 complex are taken into account. Noticeable 
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Fig. 4. Comparison of K-shell fluorescence yields, a for Fe ions 
10 < N <25. Filled circles; hfr value obtained for the first K- vacancy 
level, 10 < N < 17 (Paper III), 18 < N < 25 (this work). Filled 
squares: autostructure value obtained for the first K-vacancy level. 
Paper III. Filled upright triangle: hfr value obtained by averaging on 
all the K-vacancy levels, 10 < N < 17 (Paper III), IS < N < 25 
(this work). Filled inverted triangle: autostructure value obtained by 
averaging on all the K-vacancy levels. Paper III. Circles: Jacobs & 
Rozsnyai (1986). Squares: Kaastra & Mewe (1993). Filled diamonds: 
value recommended by Hubbell et al. (1994). Upright triangle: Sole 
(1992). Inverted triangle: Pious et al. (1992). Diamond: Bhan et al. 
(1981). Dotted circle: Singh et al. (1990). 


Fig. 3. Comparison of Auger rates for Fe ions with 18 < N < 25: 
(a) KLM/KLL and (b) KMM/KLL ratios. Filled circles: hfr, this 
work. Upright filled triangles: Chen et al. (1979). Circles: Bhalla & 
RamsdaJe (1970). Squares: Kovalik et al. (1987). Upright triangles: 
Mehlhorm & Albridge (1963). 


differences in the ratios from these two cases may be appre- 
ciated indicating sensitivity to the excitation mechanism. The 
HFR2 ratio for Fe ii is closer to the Dirac-Fock value of 
Scofield (1974), who also considered all the K-vacancy states, 
and to the solid-state experiments (Holzeret al. 1997; McGrary 
et al. 1971; Salem & Wimmer 1970; Williams 1933). 

In Figure 2, theoretical and experimental K0/K a intensity 
ratios are plotted as a function of electron number. Most exper- 
imental ratios (Rao et al. 1986; Berenyi et al. 1978; Salem et al. 
1972; Slivinsky & Ebert 1972; Hansen et al. 1970) have been 
scaled down by 8.8% in order to extract the radiative-Auger 
(5%) and 1% (3.8%) satellite contributions from the KJ3 UTA 
which theory does not include (Verma 2000). The value quoted 
in Perujo et al. (1987) does not take into account the blend with 
the K )(?5 satellite and has been corrected accordingly. With the 
exception of the values by Kaastra & Mewe (1993), theory pre- 
dicts a decrease of the ratio with N , and the theoretical scatter 
is comparable with that among the solid-state experiments of 
just under 20%. The computed results by hfr, autostructure 
and Jacobs & Rozsnyai (1986) for N = 25 are in good agree- 
ment with the bulk of the experimental values (Perujo et al. 
1987; Rao et al. 1986; Ber6nyi et al. 1978; Salem et al. 1972; 
Slivinsky & Ebert 1972; Hansen et al. 1970). Holzer et al. 
(1997) mention a possible systematic deviation in one of their 
corrections to explain the discrepancy of their data with other 
measurements. On the other hand, the Dirac-Fock ratios by 


Scofield (1974) and Jankowski & Polasik (1989) at N = 25 ap- 
pear significantly higher. We have therefore proceeded to verify 
these values by using the same code (mcdf-sal) as in Jankowski 
& Polasik (1989), and as shown in Figure 2, they are accurately 
reproduced; moreover, a similar decrease in the ratio with N is 
also obtained with this method. The spread of the different data 
sets in this comparison indicates a probable accuracy of our hfr 
transition probabilities of - 10%. 

A decrease with N is also predicted by hfr for the Auger 
KLM/KLL and KMM/KLL ratios (see Figure 3). For Fe ii, 
present results are in excellent accord with the measurements of 
Kovalik et al. (1987) and other theory (Chen et al. 1979; Bhalla 
& Ramsdale 1970) but somewhat higher than the older experi- 
mental estimate of Mehlhorm & Albridge (1963). The total K- 
shell fluorescence yields, a) K , forFe ions with 10 < N < 25 are 
presented in Figure 4. In both hfr and autostructure data sets, 
the fluorescence yields have been computed for fine-structure 
K-vacancy levels: 


= 


A r {i) 

A r (i) + A a (i) * 


( 10 ) 


In order to display the effect of the level population, we report 
in Figure 4 the value for the first K-vacancy level and the yield 
averaged over all the K-vacancy levels. As expected, the flu- 
orescence yield is independent of the population mechanism 
in the third row ions where the K /?, KLM and KMM decay 
channels become spectators. Given the complexity of level-to- 
level calculations of Auger rates, it has not been possible to 
compute fluorescence yields with autostructure for N > \1. 
The most recent experimental measurements in the solid (Sole 
1992; Pious etal. 1992; Singh et al. 1990; Bhanetal. 1981) are 
plotted along with the value recommended by Hubbell et al. 
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Fig. 5. Simulations carried out near 6.4 keV with the xstar modeling code for several values of the ionization parameter £. 


(1994). They show a scatter within 10% and are slightly lower 
than our hfr value for N = 5, i.e. 5 to 10% with respect to 
Hubbell et al. (1994) and Bhan et al. (1981) respectively. This 
is not a significant discord considering the expected accuracy 
of our fluorescence yields (~ 20%) as discussed in Paper in. 
As mentioned in the previous section, the missing participator 
Auger channels will not contribute significantly to the Auger 
widths (less than 1 %). Concerning the shake-up channels, they 
will not affect the Auger widths because they are due to mix- 
ing between valence shell configurations which have the same 
Auger widths. We have verified that this is the case in a sim- 
ple calculation ([ls](3d + 4s) 8 -> [2p] 2 (3d + 4s) 8 +£s). Data by 
Kaastra & Mewe (1993) have been corrected so as to reproduce 
the best solid-state yields compiled in Bambynek et al. (1972) 
and are therefore somewhat lower than our values. Results by 
Jacobs & Rozsnyai (1986) for N > 12 predict a substantially 
lower yield. 

In order to evaluate the impact of the new atomic data on 
Fe K-line modeling, the data sets generated in the present work 
and those in Papers I— III have been included in the xstar mod- 
eling code (Kallman & Bautista 2001). Runs for different ion- 
ization parameters, £, are shown in Figure 5. The other param- 


eters have been assigned the following values: cosmic abun- 
dances; a gas column density of 10 16 cm’ 2 ; a gas density of 
10 12 cm’ 3 ; an X-ray source luminosity of 10 38 erg/s; a power- 
law index of 1 ; and energy bins (or channel widths) of 1 eV. In 
Figure 5, one clearly sees that the UTA centroid is redshifted 
when £ goes from 0.001 to 1 and then blueshifted for f > 1 . The 
shape of the UTA changes also considerably for f > 1 where 
the ionization balance favors the first and second row ions. 

Electronic Tables 3-4 list radiative and Auger widths for 
295 energy levels; and wavelengths, A-values and fluorescence 
yields for 396 transitions. 


5. Summary and conclusions 

Following the findings of our previous study on the second-row 
iron ions (Paper III), the hfr package (Cowan 1981) has been 
selected to compute level energies, wavelengths, decay rates 
and fluorescence yields for the K- vacancy states in Fe ii— F e ix. 
The calculations have focused states populated by ground- 
state photoionization. Due to the complexity of the level-to- 
level Auger calculations, we have employed a compact formula 
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(Palmeri et al. 2001) to compute Auger widths from the hfr ra- 
dial integrals. 

The hfr centroid wavelengths for the K a and Kj3 UTAs 
in Fe 11 reproduce the solid-state measurements (Holzer et al. 
1997) to better than 1 mA. Moreover, the red shift predicted by 
hfr for the Kor lines in species with higher ionization stage is in 
accord with the EBIT wavelengths in Fe x (Decaux et al. 1995) 
thus contradicting the previous trend specified by Kaastra & 
Mewe (1993). We have also carried out extensive comparisons 
of different hfr branching ratios, namely Kor 2 /Kari, KfifKa, 
KLM/KLL, KMM/KLL and with other theoretical and ex- 
perimental data. The present ratios for Fe ii are in good agree- 
ment with the solid-state measurements, and the Ka^Kari ratio 
has been found to be sensitive to the excitation mechanism. It 
has been possible from these comparisons to estimate an accu- 
racy of -10% for the hfr transition probabilities. 

The new atomic data sets for the whole Fe isonuclear se- 
quence that have emerged from the present project are provid- 
ing a more reliable platform for the modeling of Fe K lines. 
Preliminary simulations of the emissivity of a photoionized gas 
with xstar (Kallman & Bautista 2001) have shown K a line pro- 
files and wavelength shifts that are sensitive to the ionization 
level of the gas. Further work will therefore be concerned with 
improving diagnostic capabilities by means of accurate K-shell 
photoionization and electron impact excitation cross sections. 
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